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Introduction 

In the approaching era of career astronauts and 
space workers who will man the Space Station, 
knowledge of galactic cosmic-ray interaction and 
transport in bulk matter is required to accurately an- 
alyze requirements for shielding from space radiation 
and to adequately assess the resulting radiobiologi- 
cal damage to the astronauts themselves. In previ- 
ous work (refs. 1 and 2 and references cited therein), 
a simple theory of projectile heavy-ion fragmenta- 
tion based upon an abrasion-ablation model was pre- 
sented. For the abrasion step, a quantum mechanical 
formalism based upon an optical model potential ap- 
proximation was used (ref. 3). The charge dispersions 
of the excited projectile prefragments were calculated 
by both a hypergeometric distribution and a method 
based upon the zero-point oscillations of the giant 
dipole resonance (ref. 4). The excitation energies of 
the projectile prefragments were estimated from the 
geometric “clean-cut” model (ref. 5). In the ablation 
stage of the fragmentation process, compound nu- 
cleus decay probabilities for the various particle emis- 
sion channels of the excited prefragment were calcu- 
lated by the Monte Carlo program EVAP-4, which 
was developed at Oak Ridge National Laboratory 
(ref. 6). 

In this work, the simple theory of projectile heavy- 
ion fragmentation is extended to include the effects 
of Coulomb dissociation (ref. 7) and frictional spec- 
tator interactions (FSI) (ref. 8) and then applied to 
the study of fragmentation of target nuclei by rel- 
ativistic protons and heavy ions. Target nucleus 
fragmentation is of interest to the galactic heavy- 
ion shielding problem because of its relevance to 
radiobiological damage studies and because it rep- 
resents a substantial source of experimental data 
suitable for validating our heavy-ion fragmentation 
theory. 


involved in predicting target nucleus fragmentation 
cross sections is given in figure 1. 

Abrasion Cross Sections 

Interchanging target and projectile nuclei in the 
expressions presented in reference 1 gives the cross 
section for abrading m target nucleons as 

a t \ 

m J 

where the probability of not removing a target nu- 
cleon in the collision is (ref. 1) 

P( b) = exp [-Ap a(e) /(b)] (2) 


J d 2 b [1 - P(b)] m P{b) AF (1) 



with 
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and prefragment mass number 

Ap = Ap — m 


[1 - 0(y)l 

( 3 ) 

( 4 ) 


(Definitions of symbols and abbreviations used in 
this paper appear after the references.) Methods 
for determining the appropriate projectile and tar- 
get nuclear distributions pi and constituent-averaged 
nucleon-nucleon cross sections a(e) are described in 
references 9 and 10. Values for the nucleon-nucleon 
scattering slope parameter are obtained from the ap- 
propriate parameterization for diffractive scattering, 
which is given in reference 11. As before (refs. 1 and 
2), the Pauli correlation is 


Target Fragmentation Theory 

Since the underlying physical concepts applicable 
to the abrasion-ablation fragmentation model have 
been discussed elsewhere (refs. 1 and 2), they are not 
repeated herein. Instead, a brief outline of the cal- 
culational methods, as they are used to predict tar- 
get fragmentation cross sections, is described in the 
following subsections. For more detailed discussions, 
the reader is urged to consult the references. The va- 
lidity of using the projectile fragmentation theory to 
also describe target fragmentation rests upon the un- 
derlying assumptions of the abrasion-ablation model 
(ref. 8) and the symmetry between projectile and 
target ions inherent in our optical model formalism 
(ref. 2). A flow diagram of the calculational steps 


C(y) = 0.25 exp ^ —p ^ (5) 

where 

kp — 1.36 fm -1 

Prefragment Charge Distributions 

Since the abraded nucleons consist of protons 
and neutrons, which are not identical, a prescription 
for calculating the charge dispersions of the prefrag- 
ments is needed to calculate final, isotope, and/or 
elemental production cross sections caused by the 
fragmentation process. Two such methods are used 
in the fragmentation theory described in this work. 



The first method (ref. 1) treats the neutron and pro- 
ton distributions as completely uncorrelated. The 
cross section for forming a particular prefragment of 
mass Aj and charge Zj is then given in terms of the 
cross section for abrading m nucleons a m as 



i N ) ( z ) 

&abr{Zj, Aj) — ( At ) ° m 

(6) 

where z of the original Z target nucleus protons 
are abraded along with n of the original N target 
neutrons. Note that 


Ap — N + Z 

(7) 

and 

m — n + z 

(8) 

with 

Zj — Z — z 

(9) 

and 

& 

ii 

£ 

i 

3 

(10) 


The hypergeometric distribution is based on the as- 
sumption that there is no correlation at all between 
neutron and proton distributions. Therefore, un- 
physical results such as abrading all neutrons or pro- 
tons from a nucleus while leaving the remaining frag- 
ment intact could occur. 

As an alternative to the hypergeometric distribu- 
tion, Morrissey et al. (ref. 4) proposed a charge dis- 
persion model based upon the zero-point vibrations 
of the giant dipole resonance of the target nucleus. 
In this model, equation (6) becomes 


allowed values of Zj yields unity for the dispersion 
probabilities. 

Prefragment Excitation Energies 

The excitation energies of the target prefragments 
after abrasion of m nucleons are assumed to be given 
by 

#exc = -Esurf + %SI ( 14 ) 

where i£ sur f is the surface energy calculated from 
the clean-cut abrasion formalism, and £psi> the 
contribution from frictional spectator interactions, 
is estimated by the methods in reference 8. For 
completeness, the calculational details are included 
in this work. 

From the clean-cut abrasion formalism (ref. 1), 
the surface energy term is 

E snd = A • E s (15) 

where the liquid drop nuclear model surface energy 
coefficient E s is 0.95 MeV/fm 2 , and the excess sur- 
face area of the misshapen nucleus (following abra- 
sion) is 

A = 4tt/$[1 + P - (1 - F) 2 / 3 ] (16) 

The expressions for P and F vary for different types 
of collisions (peripheral or central) and for different 
relative sizes of the colliding nuclei. The target 
fragmentation expressions for P and F are obtained 
from equations (15) through (27) of reference 1 by 
interchanging Rp and Rp. For peripheral collisions 
(Rp — Rp <b< Rp + Rp) where Rp > Rp we have 


0 _,/9 

°abr i Z v A j) = N j ( 27ra |) 1/2 ex P 1 ^2 

where the variance (dispersion) is 


Vm 

(ii) 


-= 26 K ( i ) ,/ 2 fe )(^) a+ “ r3/4 (12) 


with 

(13) 

In the droplet model of the nucleus, the coefficients 
J and Q have the nominal values of 25.76 MeV 
and 11.9 MeV, respectively (ref. 4). The rate of 
change of the number of nucleons removed as a 
function of impact parameter (dm/db) is calculated 
numerically by using the geometric abrasion model 
of reference 12. The normalization factor Nj insures 
that for a given value of Ay, the discrete sum over all 


P = o.iW/>(I- 2 )(i^) 2 

-0.12S[0.5(pv}*/2 (^ — 2 ) + 1] (~~y (17) 

and 

F = 0.75(1 — j /) 1 / 2 ( ^7^ ) 2 - <U25[3(1 - 1/) x / 2 - 1] 

"(18) 


with 


Rp 

Rp -j- Rp 


(19) 


and 


^ Rp 4- Rp 

= 1 _1 == ^P 

^ v Rp 


( 20 ) 


( 21 ) 


If the collision is central, then the target vol- 
ume completely overlaps the projectile volume 
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(b < Rp — Rp), and all the target nucleons are 
abraded. In this case, equations (17) and (18) are 
replaced by 


p = -1 

(22) 

F — 1 

(23) 


and there is no ablation of the target, since it was 
destroyed by the abrasion. 

For the case where Rp > Rp and the collision is 
peripheral, equations (17) and (18) become 


p - al26 <'“'>‘ /s ( J - ■ 2 ) (^r) ~ (;) 1/a 

*( H - 

*(¥) 3 

and 

F = 0 . 75(1 - u) 1 / 2 - 0 . 125 | 3(1 ~^ )1/2 ^ 

^ [1 ~ U - m 2 ) 3/2 ][i - (i - m ) 2 ] 1/2 I ^i-^ 3 (25) 



where the impact parameter is restricted such that 


Rp — Rp < b < Rp + Rp (26) 


For a central collision (6 < Rp — Rp) with Rp > Rp, 
equations (24) and (25) become 




1(1 _ ^ 2 ) 1/2 _ ! 



(27) 


and 


F = [ 1 - (1 - m 2 ) 3/2 J 


-0“ 


I 1/2 


(28) 


In the abrasion process, some of the abraded nu- 
cleons are scattered into, rather than away from, the 
prefragment. This scattering results in the deposi- 
tion of additional excitation energy by these so-called 
“frictional spectator interactions.” From reference 8, 
the excitation energy due to p frictional spectator 
interactions is 


&FS1 -P< E dep > (29) 


0 to m, and < E dep > is the average energy de- 
posited by any of the p nucleons. From reference 8, 
the rate that energy E (in megaelectronvolts, MeV) 
is deposited by an interacting nucleon traversing a 
distance x into the prefragment is 


dE _ _E_ 
dx 4A 


(30) 


where the mean free path is 


A = ( P<tnn ) 1 (31) 


with an assumed nuclear density of p = 0.17 fm~ 3 
and an approximate nucleon-nucleon cross section (in 
femtometers squared, fm 2 ) given by 


300 

°NN « “g- 


(32) 


Since the maximum deposited energy occurs when 
the abraded nucleon traverses a target nucleus diam- 
eter, then 


Emax — %Rt 


dE 

dx 


(33) 


where the uniform density target nucleus radius is 


Rt = r 0 A J, /3 


(34) 


and rg = 1.2 fm. In terms of the maximum deposited 
energy, the average deposited energy, obtained from 


< Edep ^ 


7T 


-1 



^dep(^) 


(35) 


is 

^ -Etfep ^ ^ g'Emax (36) 

where 0 is the angle describing the orientation of 
the trajectory of the abraded nucleon with respect 
to its diameter within the prefragment. Substituting 
equations (30) through (34) into equation (36) yields 

< E dep > « 10.2 A^ 3 (37) 

If one assumes that each abraded nucleon, on the 
average, has a 50-percent probability of being scat- 
tered back into the target prefragment, then each 
prefragment with mass number Ay and charge num- 
ber Zj which has experienced p frictional spectator 
interactions has an abrasion cross section given by 

( m ) 

a abr ( Zj j Ay , p ) = ^abr ( Zj , Ay ) (38) 


where p, the number of abraded nucleons scattered 
into the prefragment, takes on integer values from 
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with cr a br(^?> Aj) given by either equation (6) or 
equation (11). 

Ablation Factors 

Depending upon the excitation energy, the decay 
of the excited prefragment may occur by emission of 
one or more nucleons (protons or neutrons), compos- 
ites (deuterons, tritons, 3 He, or alpha particles), or 
gamma rays. The probability a^j for formation of a 
particular final fragment of type i as a result of the 
deexcitation of a prefragment of type j is obtained 
from the EVA-3 computer code (ref. 12). The final 
nuclear fragmentation cross section for production of 
the type-i isotope is then given by 

m 

flnuc ( Zj , Aj) ~ ^ ^ ^ ^ Ot-ij , Ay , p) (39) 

3 

where cr a b r (-2 J ’, Ay, p) is obtained from equation (38). 
To compare these predictions with experimental 
data, one must also include contributions due to 
Coulomb dissociation (ref. 7) to yield a cross sec- 
tion for production of a final fragment with charge 
number Z{ and mass number as 

Ai) = ^nuc [Zfa A z *) + (40) 

Elemental production cross sections are obtained by 
summing over all isotope contributions as 

a z( Z i) = ^2 GF{ Z i, A i) ( 41 ) 

A 

Results 

Cross sections for the production of carbon, 
boron, beryllium, and lithium isotopes in the frag- 
menting of carbon target nuclei by relativistic car- 
bon, neon, and iron projectile nuclei are listed in 
tables I to III. Displayed are theoretical predictions 
using both the hyper geometric (eq. (6)) and giant 
dipole resonance (eq. (11)) charge dispersion mod- 
els. The experimental results for carbon-carbon col- 
lisions (in table I) are taken from references 13 and 
14. The reference 13 cross sections are actually for 
fragmentation of the carbon projectile, rather than 
the carbon target, but are useful for comparison pur- 
poses. In general, as shown in table I, the giant 
dipole resonance model yields better agreement with 
experimental data. Tables II and III show reasonable 
agreement between theory and experimental values 
(ref. 15) for the 11 C production cross sections, but 
the lack of data on experimental cross sections for 


other isotopes precludes further meaningful compar- 
ison. For all three collisions, the theoretical cross 
sections were calculated by using the steps displayed 
in figure 1. 

In table IV, representative cross sections for the 
fragmentation of copper targets by high-energy pro- 
tons are displayed. Giant dipole resonance charge 
dispersions were not calculated, since the correct 
value to use for the proton radius is not clear. Choos- 
ing a value 1 fm) based upon the proton charge 
radius yielded unphysical values for a z (eq. (12)). In 
addition, negative values for E sur f (eq. (15)) were 
obtained with protons because of the approximate 
nature of equations (17) and (18). Therefore, we as- 
sumed that E sur f = 0 for protons and used 

£5exc — #FSI (42) 

in the p + 64 Cu calculations. Much of the dis- 
agreement between theory and experimental values 
(ref. 16), however, results from the “cumulative” na- 
ture of the experimental cross sections. These cross 
sections were obtained by counting the irradiated 
target samples after periods of time ranging from 
several hours to several days. Hence, most of the 
short-lived isotopes produced in the fragmentations 
decayed prior to counting. A more realistic method 
of comparison would be to correct the theoretical re- 
sults for these decays and predict “cumulative” cross 
sections. This comparison was not done, however, 
because cumulative cross sections are of very limited 
usefulness in radiobiological applications of the frag- 
mentation theory. 

As suggested by these results, future improve- 
ments to the fragmentation theory will center on de- 
veloping correct methods for estimating the prefrag- 
ment excitation energy. One possibility which will 
be explored is the impulsive excitation model of ref- 
erence 17. 

Concluding Remarks 

A simple abrasion-ablation model for describing 
heavy-ion fragmentation has been extended to in- 
clude frictional spectator interaction effects and ap- 
plied to target nucleus fragmentation processes. The 
model was shown to yield reasonable agreement with 
experimental isotope production cross sections ob- 
tained at Lawrence Berkeley Laboratory in colli- 
sions of relativistic carbon, neon, and iron projec- 
tiles with carbon targets. Predictions of target frag- 
mentation cross sections for proton-copper collisions 
were also made. For nucleus-nucleus fragmentations, 
the giant dipole resonance charge dispersions ap- 
peared to yield better agreement between theory and 
experiment than the hypergeometric model. The 


4 



major shortcoming of the fragmentation model, at 
this time, appears to be the methods used to estimate 
the prefragment excitation energies. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
March 27, 1986 
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Symbols 

Ap prefragment nuclear mass number 

Ai nuclear mass number for ith nucleus 

Ap projectile nuclear mass number 

Ap target nuclear mass number 

B(e ) average slope parameter of nucleon- 

nucleon scattering amplitude, fm 2 

b target impact parameter, fm 

C(y) Pauli correlation function, defined 

in equation (5) 

^dep energy deposited by nucleon under- 

going FSI 

Eex c target prefragment excitation 

energy, MeV 

Epsi frictional spectator interaction 

excitation energy, MeV 

E max maximum deposited energy, MeV 

E s nuclear surface energy coefficient, 

MeV/fm 2 

E SVlT f nuclear surface excitation energy, 

MeV 

e two-nucleon kinetic energy in their 

center of mass frame, GeV 

F function defined in equations (18), 

(23), (25), and (28) 

FSI frictional spectator interaction 

/(b) defined in equation (3) 

J droplet model coefficient 

(25.76 MeV) 

m number of abraded nucleons 

N total number of nuclear neutrons 

Nj renormalization coefficient for jth 

prefragment 

n number of abraded neutrons 

P function defined in equations (17), 

(22), (24), and (27) 

P(b) probability of not removing a single 

nucleon by abrasion 

p number of abraded nucleons con- 

tributing to FSI 


Q 

droplet model coefficient 
(11.9 MeV) 

Rp 

uniform nuclear radius of projectile, 
fm 

Rp 

uniform nuclear radius of target, fm 

ro 

uniform nuclear radius parameter 
(1.2 fm) 

u 

defined in equation (13) 

X 

distance into prefragment, fm 

y 

two-nucleon relative position, fm 

z 

total number of nuclear protons 

Zi 

total number of nuclear protons for 
ith nucleus 

z 

number of abraded protons 

z 0 

longitudinal position of target 
center of mass in projectile rest 
frame, fm 

( A ) 

\mJ 

binomial coefficient 

a ij 

probability of formation of type-z 
fragment as a result of deexcitation 
of type-/ prefragment 

az 

defined in equation (12) 

P 

defined in equation (20) 

A 

excess nuclear surface area, fm 2 

0 

orientation angle, rad 

X 

mean free path, fm 

P 

defined in equation (21) 

V 

defined in equation (19) 

tP 

collection of constituent relative 
coordinates for projectile, fm 

P 

nuclear density, fm” 3 

PP 

projectile nuclear density, fm -3 

PT 

target nuclear density, fm” 3 

a(e) 

average nucleon-nucleon total cross 
section, mb 

Vabri^j, Aj) 

cross section for production of 
nucleus of type (Z, A) by abrasion, 
mb 

&abr(^i P ) 

abrasion cross section for nucleus 
of type (Z, >1) which has had p FSI 
interactions, mb 
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<?EM 

Coulomb dissociation cross section, 


mb 


fragmentation cross section, mb 

Om 

cross section for abrading m nucle- 


ons, mb 


&NN 

defined in equation (32) 

^nuc 

final nuclear fragmentation cross 


section, mb 

°Z 

element production cross section, 


mb 
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TABLE L ISOTOPE PRODUCTION CROSS SECTIONS FROM TARGET 
FRAGMENTATION FOR THE REACTION 12 C + 12 C ->A Z + X 

[Incident kinetic energy is 2.1 GeV /nucleon} 


■ 

Isotope production cross sections, mb 

Isotope 

produced 

Hypergeometric 

model 

Giant dipole 
resonance model 

Experiment 
(ref. 13) 

Uq 

56.4 

56.4 

' 

°60.9 ± 0.6 
46.5 ±2.3 

io c 

28,9 

13.7 

4.11 ±0.22 

9 C 

8.2 

1.6 

0.539 ±0.066 

n B 

56.4 

56.4 

53.8 ± 2.7 

10 B 

23.4 

31.0 

35.1 ±3.4 

9 B 

52.3 

61.8 


8 b 

18.3 

13.4 

1.72 ±0.13 

7 b 

6.8 

1.8 


^Be 

29.2 

14.0 

5.81 ±0.29 

9 Be 

51.1 

59.9 

10.6 ±0.5 

7 Be 

46.8 

54.1 

18.6 ± 0.9 

y Li 

8.2 " : 

1.6 

0.851 ±0.082 ' 

8 Li 

18.3 

13.4 

2.18 ±0.15 

7 Li 

32.0 

55.9 

21.5 ±1.1 

6 Li 

69.3 

86.1 

30.0 ±2.4 


ft Value from reference 14. 




TABLE II. ISOTOPE PRODUCTION CROSS SECTIONS FROM TARGET 
FRAGMENTATION FOR THE REACTION 20 Ne + 12 C A z + X 

[Incident kinetic energy is 1.05 GeV/nucleon] 


Isotope 

produced 

Isotope production cross sections, mb 

Hypergeometric 

model 

Giant dipole 
resonance model 

Experiment 
(ref. 15) 

1T C 

65.3 

65.3 

80.3 ±1.4 

10 c 

32.7 

16.6 


9 C 

9.4 

2.0 


n B 

65.4 

65.4 


i°B 

27.3 

35.4 


9 B 

60.2 

70.6 


8 b 

20.7 

15.6 


7 b 

8.0 

2.2 


roge 

32.8 

16.7 


9 Be 

57.3 

67.8 


8 Be 

54.1 

60.0 


9 Li 

9.4 

2.0 


8 Li 

20.7 

15.6 


7 Li 

57.4 

64.5 


6 Li 

79.7 

99.0 
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TABLE IV. REPRESENTATIVE ISOTOPE PRODUCTION CROSS SECTIONS FROM 
COPPER FRAGMENTATION IN THE REACTION p + 64 Cu -> A z + X 

[Incident proton kinetic energy is 3,9 GeV] 


Isotope 

Isotope production cross sections, mb 

produced 

Hypergeometric model 

Experiment (ref. 16) 

B3 Cu 1 

66.2 

14.3 ±0.7 

62 Cu 

11.8 


61 Cu 

14.8 

11.5 ±0.4 

^Cu 

13.1 

1.92 ± 0.6 

83 Ni 

54.8 


62 Ni 

20.0 


81 Ni 

23.8 


eo N i 

69.9 


59 N j 

33.4 


58 N i 

10.4 


62^ 

8.0 

0.51 ±0.04 

61 Co 

7.5 

4 ±0.08 

60 Co 

14.9 

8.34 ±0.09 

59 Co ! 

22.2 


58 Co 

29.3 

20.6 ± 0.2 

57 Co 

23.0 

16.9 ±0.1 

56 Co 

15.6 ! 

5.6 ±0.3 

5g __ 1 

3.9 ; 

1.33 ± 0.24 

™Fe 

1.9 ! 

1.83 ± 0.22 

^ 8 Mn 

6.7 | 

2.58 ±0.08 

54 Mn 

15.1 i 

13.8 ±0.2 

52 Mn 

7.6 

6.26 ±0.08 
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Figure 1. Flow diagram of steps in target fragmentation calculation. 
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